Enhancement of Critical Temperature of a Striped Holographic Superconductor 
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We study the interplay between the stripe order and the superconducting order in a strongly 
coupled striped superconductor using gauge/gravity duality. In particular, we study the effects ol 
inhomogeneity introduced by the stripe order on the superconducting transition temperature beyond 
the mean field level by including the effects of backreaction onto the spacetime geometry in the dual 
gravitational picture. We find that inhomogeneity enhances the critical temperature relative to its 
value for the uniform system. 
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I. INTRODUCTION 

One of the differences between conventional supercon- 
ductors and high temperature superconductors is that 
the normal states of the conventional superconductors 
are well described by Fermi liquid, whose only (weak 
coupling) instability is to superconductivity. By con- 
trast, the normal states high temperature superconduc- 
tors, such as cuprates and iron pnictides, are highly cor- 
related and thus, exhibit other low temperature orders 
which interact strongly with superconductivity. One of 
the prominent orders is the unidirectional charge density 
wave "stripe" order PHD- 

It is therefore important to understand the nature of 
the interplay between superconductivity and the stripe 
order in the presence of strong correlation. In this article, 
we study the effects of the stripe order on the supercon- 
ducting transition temperature of the strongly correlated 
superconductor using holography or gauge/gravity cor- 
respondence. In gauge/gravity duality, the strongly cou- 
pled condensed matter systems are mapped to a weakly 
coupled Einstein-Maxwell-scalar theory on black hole 
spacetimes with negative cosmological constant, or the 
so-called anti de-Sitter (AdS) black holes. Just like the 
normal states of high temperature superconductors, AdS 
black holes feature numerous types of instability that 
lead to the formation of scalar 'hair ' @, HI (which cor- 
responds to superconductivity), striped phases and 
nematic phases Ultimately, we would like to study 
the system where both the superconducting order and 
the stripe order emerge dynamically, however, since our 
focus in this article is the effects of the stripe order on 
the critical temperature, we will follow Refs. [ll| and [TU 
where the inhomogeneity is introduced via a modulated 
chemical potential, with wavenumber Q. These articles 
studied the holographic striped superconductor by ne- 
glecting the backreaction of the electromagnetic field on 
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the spacetime geometry, which in field theory language 
corresponds to the mean field treatment of the system. 
To improve upon their results, we would like to consider 
the effects of fluctuations on the strongly coupled striped 
superconductor by studying the backreacted spacetimes. 



Our result shares similar qualitative features as the re- 
sult of Ref . [IH , which study the effects of local inhomo- 
geneity on the critical temperature within the framework 
of BCS theory. This similarity is of course not surpris- 
ing, as one can reproduce the weak coupling result from 
the strong coupling result by simply taking the limit in 
which the scaling dimension of the superconducting or- 
der becomes the scaling dimension of a free scalar field 



We find that the critical temperature for the forma- 
tion of the scalar hair is maximum in the case of AdS 
Schwarzschild black hole, which is the solution to the 
Einstein equation when backreaction is neglected pH[I3 [. 
This means that the superconducting transition tempera- 
ture obtained by including fluctuations is lower than that 
obtained in the case where fluctuations are neglected. In 
the field theory side, this can also be understood using 
a Ginzburg-Landau type argument: introducing fluctua- 
tions costs us free energy and thus lowering the critical 
temperature. Usually, the corrections due to fluctuations 
are small enough that the qualitative behavior of T c as a 
function of other parameters remains the same. However, 
when the system is inhomogeneous due to the presence 
of stripe order, the effects of fluctuations at the regime 
where the inverse fluctuation length scale is smaller than 
the wavenumber Q are rather drastic. Starting at Q = 0, 
when we increase Q, we see that T c exhibits a steep 
jump as we turn on the modulation and after reaching a 
maximum, it monotonously decreases as Q increases and 
asymptotes to a constant value at Q — > oo. As the crit- 
ical temperature at finite Q is larger than the values at 
Q = and Q — > oo, which correspond to the homogenous 
limits, we find that the critical temperature is enhanced 
by the presence of stripe order. 
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II. SET-UP 



where 



We are interested in studying a strongly coupled 
striped superconductor using the gauge/gravity dual- 
ity. To this end, consider a U(l) gauge potential A a 
and a scalar field ip charged under this potential liv- 
ing in a spacctimc with negative cosmological constant 
A = — 3/L 2 . The scalar field is dual to the scalar order 
parameter of the superconductor, i.e., the condensate, 
while the U(l) gauge field is dual to the four-current in 
the condensed matter system. For simplicity, we shall 
adopt units in which L = 1, 16itG = 1. 

To study the strong coupling regime of the supercon- 
ductor, we only need to study the gravity theory at the 
classical level. In particular, we are interested in find- 
ing solutions to the classical equations of motion whose 
boundary values are related to the parameters of the su- 
perconductor. 

The action for this system is 



R + &-\f 2 -\D a if -m 2 \^ 



(1) 

where, D a = d a - iqA a , F ab = d a A b - d b A a and a, b G 
{t, r, x,y}. Here 

The field equations consist of the Einstein equations, 



Rab 



9abR — 3g a fc = X^ob , 



(2) 



where the stress-energy tensor is 



F ac F b ' 



1 



-^9abF Cd F cd 



Da^D^Y + (a O b) - g ab [\D a ip\ 2 + m 2 \^< 



the Maxwell equations 
1 



where the U(l) current is 

J a = -#*£> a V - c.c] , 
and the Klein-Gordon equation for the scalar field, 

1 - - ' ab 



-9 



D a (V=gg aD D b il;)+m^ = 0. 



(3) 



(4) 



(5) 



(6) 



III. ABOVE THE CRITICAL TEMPERATURE 

Above the critical temperature T c , the scalar field van- 
ishes (tp = 0) and thus, the Einstein-Maxwell equations 
simplify to 

R\ + Z5\ = \t\, d b (V=9~F ab ) = 0, (7) 



- \6\F cd F cd . 



(8) 



We are interested in finding a static black hole solution of 
flat conformal boundary which is sourced by an external 
modulated chemical potential /it(x), where \i is a given 
spatially dependent function. Let (x, y) be the Carte- 
sian spatial coordinates of the two-dimensional confor- 
mal boundary. We concentrate on the case in which (i 
only depends on one of the coordinates, which is chosen 
to be x. Expanding in Fourier cosine modes, we write 



H(x) 



oo 

H S n cos nQx . 

n=0 



n=0 



(9) 



For definiteness, we consider the case in which only the 
first couple of modes are no n- vanishing, setting 



5 = 1-6, S 1 = 6, 6 n = (n > 2). 



(10) 



To find a solution to the Einstein-Maxwell equations (J7J , 
consider the metric ansatz 



j 2 2 - 

as = —r e 



rir 2 



p [e^dx 2 + e^dy 2 ] , 
(11) 

where a, ft, and 7 are functions of (r, x). The boundary 
is at r — > 00 and the horizon is at r = r+, where r + is 
an arbitrary parameter. For a flat conformal boundary, 
we require a, f3, 7 — > 0, as r — > 00, and in fact, we find 
a - 0(r~ 3 ) while ^ and 7 - 0(r~ 4 ). 

For the U(l) potential, we fix the gauge such that A r = 



A, 



A y = and A t = A t (r,x) with A t = at the 



horizon (for a finite norm, A a A a < 00), whereas at the 
boundary, 



A t (r,x) 



H(x) 



(12) 



We shall solve the Einstein-Maxwell equations ([7]) per- 
turbativcly in [i 2 . This expansion is valid for large black 
holes (or small chemical potential), or more precisely for 



< 



(13) 



Expanding in the small dimensionless parameter 
we have 



4 0) + I ^ 



r+ 

/' 

r 



A^ 



P = f3 (0) + 
7 = 7 (0) + 



/' 



(14) 
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and consequently, the expansions of the metric, Ricci ten- 
sor and U(l) field strength and stress-energy tensor, re- 
spectively, 



(o) . 

9ab = <?ab + 



Rab 
F a b 



9$ 



R ab+[- 



R 



(i) 



F 



(0) 



a b 



Tab= ^ T n 



-(0) 



(15) 



Explicitly, the non-vanishing components of the 
zeroth-order electromagnetic stress-energy tensor are 



j-{0) t 



.r(°)f =-— 

v 4 

_ T (0)x = l_ 
I x 4 



c2 

-^+£ 2 

h z 

f 2 

h tz 



T (0)x _}_ T (0)z _ _^_f c 

' z ~h x ~ 2h xz 



(23) 



given in terms of the components of the electric field 
c x = — Ai sin Qx , 



At zeroth order, the Einstein-Maxwell equations read 



R 



(0)a 



+ 35% = o , d b (V-g {0) F {0)ab ) =0. (16) 



The Einstein equations decouple and are solved by the 
AdS Schwarzschild black hole 



= h = l- 



P±) 3 , /}(0) =7 (0) =0 . (17) 



To solve the Maxwell equations, it is convenient to intro- 
duce the coordinate 



(18) 



so that the boundary is at z = and the horizon at z = 1 . 
Expanding the U(l) potential in Fourier modes 



A[ ^ = \i S n A n (z) cosnQa 



we deduce the mode equations 



r\h{z) 



A n (z) = , 



(19) 



(20) 



to be solved together with the boundary conditions 



A n {0) = 1, A n {\) = 0. Here, h(z) = 1 



(Eq. (HB) 



and ' denotes a derivative with respect to z. For n = 0, 
we obtain 



Ao{z) = l-z. 



(21) 



For n > 1, a good analytic approximation to the solution 
is given by 



sinh 



A n (z) 



(1-z) 



sinh^ 



(22) 



With the choice of boundary conditions (|T()|) , the lowest- 
order stress-energy tensor has modes with n < 2 due 
to the fact that it is quadratic in the U(l) potential. The 
same should be true for the first-order corrections to the 
metric. 



£ z = (1 - S)Mq (z) + 6 A[(z) cos Qx. 
To solve the Einstein equations at first order, 



R 



(l)a 



T 



(0)a 



6 ■ 



(24) 



(25) 



we set 








= o$\z)- 










7 (i) 


= 7^)- 


h 7 P 



„(!)/ 



(26) 



We obtain five non-vanishing components for each set of 
functions {a[ X \ /3^\ % }, where i = 0, 1,2. Of the five 
equations, only three are independent and can be solved 
analytically for the three corresponding metric functions. 
After some algebra, we obtain the following system of 
equations for the modes of the metric functions. 
For the Fourier zero modes, we obtain 



.X 1 ) 



ti 



(i)' 



(^5 2 A\-h (2(1 - S) 2 A' 2 + 5 2 A 12 )) 



3 



8h 2 



(1)" 
7o 



} 2 z 2 S 2 A 2 
4r 2 + h 2 

fz 2 $ 1 A\ 
Ar 2 + h 2 



0. 
0. 



(27) 



We solve these equations by requiring that the functions 
be regular at the horizon (z = 1) and vanish sufficiently 
fast at the boundary {z = 0). We obtain 



Q 2 5 2 



dz' / dz" 
Jo 



{z"f A\ 
h 2 



2 



(28) 



(29) 



4 



a£\z) = ^fa£\z')dz', 



(30) 



where 



4%) = 2(1-*)^+^-^^-^ 



1/2 



For the Fourier first modes, we obtain 



A\ (i)/4/t- zti 




z 3 £(l - 
2/i 

- - ^ M'" 



(i)" 



Ti 



The second equation readily yields 



Q 2 



= 0. 

= o : 

0. 



M 1} (z) = 0. 



(31) 



(32) 



(33) 



The second equation yields 



4r 



We note that $p = 

Eliminating a 2 between the other two equations, we 
obtain, as before, a third-order differential equation for 
72 , from which we deduce the possible near horizon be- 
havior, 1 - z, (1 - z) ln(l - z), and (1 - z) 1+2q2 ! / 3r + . As 
before, we fix the three integration constants by demand- 
ing 7 «(0) = 0, 7^'(0) = 0, and 7 «" < 0(1/(1 - z)) 
at the horizon (z = 1). The second boundary condition, 



dz 



} (i) 



,(*') 



'\2 



dz"q± 

h 



(35) 



ti) 



at the bound- 



together with Eqs. (HU, ensure a 
ary. The third boundary condition is necessary for the 
existence of a well-defined temperature (surface gravity) , 
resulting in (1) = 0, on account of the third equation 
in (Ml). 

The equations for the various modes can be solved nu- 
merically subject to the boundary conditions outlined 
above. We have plotted o;„ (n = 0,1,2) in Fig. [T] for 
representative values of Q, whereas /% (n = 0, 2; it van- 
ishes for n = 1) is plotted in Fig. [21 and 7^' is plotted 
in Figs. 131111 andO for n — 0, 1 and 2, respectively. 



By eliminating between the other two equations, we 
obtain a third order differential equation for 7^. Then 

the possible behavior of at the horizon is found to 
be a linear combination of 1 — z, (1 — z)ln(l — z), and 
(1 — z) 1+( 3 / 6r +. We fix the three integration constants 

by demanding 7 f ) (0) = 0, 7 1 1) '(0) = 0, and -f{ 1} " < 
0(1/(1— z)) at the horizon (z = 1). The second boundary 

condition, together with Eqs. (|3"2"]l . ensure ap ~ z 3 at 
the boundary. The third boundary condition is necessary 
for the existence of a well-defined temperature (surface 

gravity), resulting in a± (1) = 0, on account of the third 
equation in (|3"2")l . 

Finally, for the Fourier second modes, we obtain 



0.2 0.4 0.6 0.8 



= 3.0 ; 5=0.5 
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ti 



(1)' 



72 



(1)// 



Q 2 (16 h 



8h 2 

6 2 Q 2 z 2 A 2 
+ Ar\h 2 

5 2 z 2 A\ 



Ar 2 h 2 



0. 
0. 

0. (34) 



FIG. 1: a ( o> (black), af 5 (red) and o$> (blue) for 
8 = 0.5, and Q/r + = 0.1, 1, 3 and 5. 



Note that, the (5^ and components of metric per- 
turbations are sourced by x-component of electric field 
(f24|) . which vanishes at both small and large Q. The 
functions ^'s and 7^ s depend on Q via two terms: 
a direct proportionality factor Q 2 and area under the 
functions A\/h or z 2 A\/h 2 . The first factor vanishes 
at Q — > 0, while the integrals vanish at Q — > 00, due 



Q 



to Ai 

1 sinh Q 

/3n and 7„ (n 



(1 — z) near the horizon. Consequently, 
0, 1, 2) are very small in both limits 



5 



Q <C r + and Q 3> r+. These functions are more sig- 
nificant in the intermediate range 2 < Q/r + < 3 and 
decay rapidly on both sides, but even when they reach 
their maximum, they remain well below unity (see Figs. 
[2][5]). Thus, their contribution to physical quantities is 
negligible in the entire range of Q. 

6=0.5 



5=0.5 



0/006 



0.0 0.2 



0.4 0.6 



FIG. 4: 7 f ) for 6 = 0.5, and Q/r+ = 0.1 (red), 1.0 
(green) and 2.0 (blue). 



5=0.5 



FIG. 2: f3 { 1} = -4 1J for 6 = 0.5, and QL 2 /r+ = 0.1 
(red), 1.0 (green), 2.0 (blue) and 8.0 (orange). 



5=0.5 



- 0.002 



FIG. 3: 7^ for S = 0.5, and Q/r+ = 0.1 (red), 1.0 
(green), 2.0 (blue) and 8.0 (orange). 

Next, we discuss the behavior of al (n = 0,1,2) 
which are physically important because they determine 
the temperature. Indeed, the Hawking temperature at 
first perturbative order is 



T 



3r+ 

47T 



1 



(36) 



Since aiP(l) = for n > 1, we have 
Q W ( l ) = a W(l) = 5o(1) 



24 



2(1 - 5) 2 +<5 2 ^1' 1 2 - 3^ 



„(!)' 



24 



(37) 



FIG. 5: 7^' for 5 = 0.5, and Q/r+ = 0.1 (red), 1.0 
(green) and 2.0 (blue). 



We can calculate these functions analytically in the two 
important limits: Q — > and Q — > oo. 

In the limit Q — > 0, we obtain the analytic expressions 



„(1) 



(1-Sf + 



6 2 



4(1 + z + z 2 ) 



O 



\ rl 



2(1 + z + z 2 ) 

S 2 z 3 
8(1 + z + z 2 ) 



(1 _ z fQ 2 /3rl + Q ^Opj 



(38) 



At Q — (or equivalently, 5 = 0), we recover the ex- 
act Reissner-Nordstrom solution representing the homo- 
geneous system 

e = e 



1 



(39) 



6 



where we used — + + from Eq. (j26|) 
with Q = 0. 

We recover the Schwarzschild solution, which is the 
solution in the probe limit, in the limit fi — > 0. As we 
increase ii, we move further away from the probe limit 
and the effects of back reaction to the metric become 
more pronounced. We reach extremality at n/r + = 2^3, 
which gives the upper bound for /i. 

It should be noted that the convergence to the homoge- 
neous system is not uniform. At the horizon, an\l) — > 
for n = 1,2, and therefore aS 1 ' docs not converge to 
its homogeneous counterpart. In other words, the limits 
Q^O and z — > 1 do not commute. It follows that there 
is a discontinuity in the temperature which depends on 
the behavior of a„ at the horizon. From Eq. (|36[) in the 
limit Q — >• 0, we obtain 



T 



3r+ 

47T 



^ 2 ((l-S) 2 + S 2 /2) 



12rj 



(40) 



which is valid for small Q. Comparing this result with the 
homogeneous case, which is recovered by setting (5 = 0, 
we obtain an enhancement in temperature upon turning 
on modulation 



AT 



Ts=o 



- 1 



i2, T H 2 -T 



(41) 



IV. THE CRITICAL TEMPERATURE 

The Klein-Gordon equation for a static scalar field 
ip(z, x) of mass m and charge q reads 

9i {V=99 U dii>) + (q 2 A 2 t -m 2 )yj = 0. (45) 

i=z,x ^ 

The mass is related to the conformal dimension A of the 
superconducting order parameter by 

™2 



A (A - 3) 



(46) 



To solve the wave equation, we expand ip m a Taylor 
series around x = 



1>{z,x) = { ^z A [F (z) + x 2 F 1 (z) + 



(47) 



The series is expected to converge rapidly at the two ends 
(Q — > and Q — > oo), where the ^-dependence of tp is 
mild. For intermediate values of Q, the x-dependence is 
more significant, but the zero mode Fq still dominates. In 
this regime, the error in the numerical analysis can be re- 
duced to the desired accuracy by including higher modes 
of the Fourier expansion. We shall leave this to future 
work, as our focus here is the temperature enhancement 
at small Q, where the calculation involves slowly varying 
functions of x. 

We obtain the leading order wave equation 



with a maximum enhancement for 6 = |. The change 
in temperature is discontinuous, but this is an artifact of 
keeping only the first order in perturbation theory. This 
change in the temperature is expected to become smooth 
(yet remain steep) as higher orders in the perturbativc 
expansion are included. 

On the other hand, in the Q>r+ regime, the contri- 
bution of Ai becomes exponentially small, and all func- 
tions except become negligible. In this regime, we 
have 



2 ^3 



4(1 + z + z 2 ) 



(42) 



So in the Q — > oo limit, we recover another exact 
Reissner-Nordstrom solution, albeit with less charge den- 
sity, 



-^i-(i + '''VV + '' 2 VV 

\ 4ri / 4rl 



(43) 



This coincides with the homogeneous solution (|39p if 5 = 
0, as expected. 

We then deduce the temperature for large Q to be 
given by 



3r+ 

47T 



H 2 (l-5f 



(44) 



K 



_A 

z 2 'h L 



2(A-1) h 

+T 

z h 



(A - 3)(/i — 1) + zh 



Ft) 



q 2 (A[ 0) ) 2 



F = 0. 



where 



= h 



' J- v 



„(!) 



.,(1) 



t(°), 



(48) 



x=0 

(49) 



and Ar(z) = Ai o) ( Z ,0). 

To solve this wave equation, it is convenient to fix the 
parameters {A, 5, q, Q/r + } and find the eigenvalue for 



(50) 



x= q ±. 

The critical temperature is then determined by Eq. 
and by the value of T/r+, which is given by (|3l)|) . Wc 
obtain 



To 
qii. 



3 

47T 



-A a <i) (1) 



A 



(51) 



At small Q, we obtain from Eq. (|41j) the enhancement in 
critical temperature 



AT C 

~t7 



A>(*-- 

I2q 2 \ 2 



(52) 
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which vanishes at the probe limit (q —> oo) and becomes 
significant away from it. However, we need to be cautious 
in taking the small q limit, as this is only a first-order 
result, which is 0(l/q 2 ). 



The wave equation is solved numerically subject to the 
boundary conditions F ~ z A at the boundary and the 
demand of regularity at the horizon (F (l) < oo). The 
results are shown in Figs. |H] and [7J for A = f and 2, 
respectively. In each case, we have chosen the other pa- 
rameters so that the curves asymptote to the same tem- 
perature as Q — > oo. We note that all curves exhibit a 
jump at Q = + , showing the enhancement of the critical 
temperature once modulation is switched on, in agree- 
ment with our analytic result (|52"T) . As Q increases, the 
critical temperature decreases monotonically. The jump 
vanishes in the probe limit which is obtained for fj, = 
(Schwarzschild black hole). For any given Q, the critical 
temperature attains its maximum value at this limit. Put 
differently, back reaction to the metric lowers the criti- 
cal temperature. Correspondingly, in the dual boundary 
system, quantum fluctuations result in a reduction in the 
critical temperature for a given modulation vector Q. 



V. SUMMARY AND OUTLOOK 

In this article, we have studied the effect of inhomo- 
gencity on the superconducting transition temperature of 
the strongly coupled striped superconductor beyond the 
mean field level, by including backreaction of the elec- 
tromagnetic field on the geometry of spacetime in the 
dual gravitational picture. We found that as we turn on 
the modulation, the critical temperature exhibits a steep 
jump. After that, as we increase Q, the critical tem- 
perature decreases until it reaches the asymptotic value. 
In other words, we found an enhancement of the critical 

A=2 



if 

0.060 r 



0.050 - 



0.045 - 



0.0 0.5 1.0 



1.5 2.0 2.5 3.0 , 



Q.22 r 
0.21 
0.20 

0.19 M 

: / 

0.18 ; 

: / 

0.17 '-> 

- / 

0.16 - 



A=l 



FIG. 7: From top to bottom: T c vs. Q for A = 2 and 
(5,q 2 ) = (0.3, oo), (0.2,7.92) and (0.1,4.22). Parameters 
are chosen so that curves asymptote to T c / (qfi) = 0.041 
as Q — > oo. 



\ temperature due to inhomogeneity that comes from the 

/ \\ stripe order. 

The discontinuous jump we see here is an artifact of 

\\ only keeping the zero mode of the scalar field in the cal- 

\\ culation and we expect that as we include the higher 

\ modes, this jump will become smooth but yet steep. It 

will be interesting to study whether the maximum of T c 

0.15 : ~ 1~ ~ Z.Z.Z.Z. corresponds to the value of Q being the inverse of super- 

F conducting correlation length scale as is seen in the BCS 



o 



result. 



FIG. 6: From top to bottom: T c vs. Q for A = 1 and 
(S,q 2 ) = (0.3, oo), (0.2,0.411) and (0.1,0.190). 
Parameters are chosen so that curves asymptote to 
T c /(qfi) = 0.149 as Q -> oo. 
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